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A sufficiently light scalar field slowly evolving in a potential can account for the dark energy that 
presently dominates the universe. This quintessence field is expected to couple directly to matter 
components, unless some symmetry of a more fundamental theory protects or suppresses it. Such a 
coupling would leave distinctive signatures in the background expansion history of the universe and 
on cosmic structure formation, particularly at galaxy cluster scales. Using semi-analytic expressions 
for the CDM halo mass function, we make predictions for halo abundance in models where the 
quintessence scalar field is coupled to cold dark matter, for a variety of quintessence potentials. We 
evaluate the linearly extrapolated density contrast at the redshift of collapse using the spherical 
collapse model and we compare this result to the corresponding prediction obtained from the non- 
linear perturbation equations in the Newtonian limit. For all the models considered in this work, if 
there is a continuous flow of energy from the quintessence scalar field to the CDM component, then 
the predicted number of CDM haloes can only lie below that of ACDM, when each model shares the 
same cosmological parameters today. In the last stage of our analysis we perform a global MCMC 
fit to data to find the best fit values for the cosmological model parameters. We find that for some 
forms of the quintessence potential, coupled dark energy models can offer a viable alternative to 
ACDM in light of the recent detections of massive high-z galaxy clusters, while other models of 
coupled quintessence predict a smaller number of massive clusters at high redshift compared to 
ACDM. 



I. INTRODUCTION 

During the last decade, a vast amount of cosmolog- 
ical data has been collected, ranging from measures of 
the Hubble constant (H n ) [1] to precision measurements 
of the Cosmic Microwave Background (CMB) [2], obser- 
vations of Type la Supernovae (SNIa) [3] and Baryon 
Acoustic Oscillations (BAO) [4]. These datasets are con- 
sistent with the standard ACDM cosmological model, 
which consists of a cosmological constant (A) plus Cold 
Dark Matter (CDM). However this model suffers from 
several short-comings, namely the fine-tuning and coin- 
cidence problems (see e.g. [5] for a review). 

These issues surrounding the cosmological constant 
have prompted investigation into alternative models [5- 
7], such as quintessence [8-10], whereby the Dark Energy 
(DE) component is dynamic and is identified with a light 
scalar field (f>(t), slowly evolving in a potential V{(j)) (for 
a comprehensive review see [11]). Since this quintessence 
field must be sufficiently light to drive an accelerated ex- 
pansion, we should expect it to couple to other matter 
forms in the universe, unless some symmetry of a more 
fundamental theory protects or suppresses it. Whilst a 
coupling to baryonic matter is tightly constrained by lo- 
cal gravity tests [12], interactions within the dark sector 
have not been ruled out by observations. A number of 
previous studies examined the cosmological implications 
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of interacting DE-CDM cosmologies (see e.g. [13-16] and 
references therein) while analysis of CMB and BAO have 
placed constraints on various coupled quintessence mod- 
els [17, 18]. These studies have shown that dark sec- 
tor interactions can leave distinctive signatures in the 
background expansion history of the universe and on cos- 
mic structure formation [19-22], at large (galaxy cluster) 
scales in particular. 

Intriguingly, a number of massive, high redshift clus- 
ters have recently been detected [24-26] which should be 
extremely rare in the standard ACDM model. It is possi- 
ble that coupled quintessence cosmologies (hereafter CQ) 
may alleviate these tensions by modifying the expansion 
history and perturbation growth [27]. It is therefore im- 
portant to understand structure formation in CQ models. 

In this paper we consider three different quintessence 
potentials and examine how the strength of the DE- 
CDM coupling (which we take to be constant) affects 
the abundance of CDM haloes. We present results for 
three commonly used analytic forms for the CDM halo 
mass function, Jenkins et al. [28], Press-Schechter [29] 
and Sheth-Tormen [30] . How the dark energy scalar field 
behaves in highly non-linear regimes has yet to be un- 
derstood and properly modelled. Indeed, the evaluation 
of the linearly extrapolated density contrast at the red- 
shift of collapse, <5*(z), in CQ cosmologies has recently 
been the subject of some debate [31]. To get a feel- 
ing for the uncertainty when modelling the non-linear 
behaviour of the scalar field, we consider two different 
methods for evaluating S+(z), (i) the spherical collapse 
model applied to CQ, and (ii) the prediction obtained 
from the non-linear perturbation equations in the New- 
tonian limit (see Ref. [31]). We first assume that all the 
CQ models share the same cosmological parameters as 
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ACDM at z = 0. Whilst this facilitates an easier com- 
parison between the models, as emphasised by Ref. [32] 
for the case of uncoupled quintessence, it is unlikely that 
they are consistent with observational data. Therefore in 
the final stage of our analysis we perform a global fit using 
the CosmoMC package [33], a Monte Carlo Markov Chain 
(MCMC) code, to find the best fit values for the cos- 
mological parameters such that each model satisfies the 
CMB+BAO+SNla+Zfo datasets. Similar analysis have 
been carried out for a large number of non-interacting 
dark energy models [34, 35]. 

In Section II wc introduce the quintessence models 
under investigation. Here we study the evolution of 
the background and perturbed quantities, pointing out 
the distinctive features of each model. We discuss the 
halo mass functions that we choose to work with in 
Section III, paying special attention to the linearly 
extrapolated density contrast at the redshift of collapse 
5* (z). In Section IV we present our predictions for halo 
abundance including a global MCMC fit for each model 
and wc conclude in Section V. 



II. COUPLED QUINTESSENCE COSMOLOGIES 

Whilst in any multi-component system the total 
stress-energy momentum tensor of the Universe must be 
conserved [36], 

i 

the stress-energy momentum tensor for each individual 
species i may in general however, not be conserved. We 
consider cosmologies comprised of radiation (7), CDM 
(c), baryons (6) and the quintessence field {<p). Neutrinos 
are assumed to be massless, forming a fraction of the 
radiation content. Then, subject to the above constraint, 
we may consider a dark sector coupling [13] 

V^)„ = Q W ,, (2) 

where the source terms Q(i) v represent the interaction. 
We adopt natural units and set 87tGn = = 1. If 

the quintessence field is coupled only to CDM, the con- 
straint Eq. (1) requires that Q(<f,) v = —Q( c ) v - Our set of 
cosmologies are described by a Lagrangian of the form 

C = - l -d»4>d^ - V{(j>) + /Imattcrfc, 4>] , (3) 

where \i are the matter fields (baryons, cold dark matter, 
neutrinos). Here, we consider a model in which only cold 
dark matter couples to the quintessence scalar field <f> and 



in which the mass of the dark matter field is a function 
of <j>: m — m(<j)). The choice m(<j)) specifies the source 
term, which as is customary, wc define as 

Q( C ) V = -/3pcd„(p , (4) 

with 

An immediate consequence of Eq. (5) is that the mass 
of the CDM matter particles no longer scale as a~ 3 and 
hence (with j3 constant) 

Pc = Pc fi a- 3 e^ . (6) 

A coupling of the form Eq. (5) can arise for instance 
after a conformal transformation of a two-metric Brans- 
Dicke theory [37] which leaves radiation and baryons un- 
coupled. Radiation remains uncoupled since its stress- 
energy momentum tensor is traceless, whilst a coupling 
to baryonic matter is tightly constrained by local grav- 
ity tests [12]. We will assume that the coupling j3 in 
Eq. (5) is constant throughout this work. The back- 
ground and linear perturbation evolution for this set of 
coupled dark energy models has been widely studied for 
constant /3, [10, 13] as well as with regard to the effects 
on structure formation [19, 21, 23]. More recently, N- 
body simulations involving time varying couplings f3(t) 
have revealed several interesting effects in the non-linear 
regime of structure formation [20]. 

The zero-component of Eq. (2) (where Q( c )o = ~PPc<j>) 
provides the conservation equations for the energy den- 
sities of each species. For Qt c )o > (< 0), the direction 
of energy transfer is from the CDM (quintessence field) 
to the quintessence field (CDM). The direction of en- 
ergy exchange then depends upon the evolution of the 
dark energy scalar field <f> and on the sign of the cou- 
pling, ft. Once a non-zero coupling of the form Eq. (5) 
is introduced, the evolution of the field is driven by the 
properties of the effective potential, V c a — V{4>) + p c {4>), 
which may possess a minimum about which the field may 
undergo damped oscillation. Hence, the scalar field ve- 
locity <j) may switch sign, reversing the direction of energy 
exchange. Such sign changeable interactions have been 
studied in the past for models where the role of dark 
energy is played by a decaying vacuum energy [38, 39]. 

A. Quintessence potentials 

In this study, we will be concerned with three different 
forms for V{(f>) which we briefly discuss below. These 
potentials have the form V = V(,f(<p; A„), where Vo is a 
constant and /(</>; A„) is a function of the field (j>, with the 
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FIG. 1. Left panel: The evolution of the quintessence field equation of state parameter = P$/ ' p$ for the three potentials for 
a variety of coupling strengths /3. Right panel: The redshift evolution of the Hubble parameter H(z) in the 2EXP and SUGRA 
models. In the lower panel we compare to the evolution of H(z) in ACDM by plotting the ratio. For the AS model, H(z) is 
identical to ACDM over the redshift range of interest, for all coupling strengths considered in this study. 



quintessence model fixed by choice of the self coupling pa- 
rameters A„ and the strength of the DE-CDM coupling 
(3. We numerically integrate the background equations of 
motion using a version of the CAMB [40] computer software 
that we have modified to include a quintessence scalar 
field coupled to CDM. To facilitate an easier comparison 
between each model, the parameters X n in the potential 
are kept fixed whilst Vq and the CDM density p c are var- 
ied (for a given j3) such that each cosmology evolves to 
the ACDM WMAP7 Maximum Likelihood (ML) cosmo- 
logical parameter values [2]. The same initial conditions 
<pi = 10~ 4 and <pi — for the field are assumed for each 
model. 

The resulting cosmology may not be consistent with 
observational datasets and so in Section IV B we per- 
form a global MCMC fit for each quintessence model 
to determine j3 and the values of other cosmological 
parameters. 

Double exponential (2EXP) - This potential is an ex- 
ample of a scaling solution and requires no fine tuning of 
the initial conditions [41]. It is given by 

V(0) = y o [e Al0 + e A ^]. (7) 

It is possible that potentials of this type could arise in 
string theory as a result of Kaluza-Klein type compact- 
ification [41]. Throughout this study, we set the param- 
eters in this potential to Ai = 20.5, A2 = 0.5 to provide 
~ -0.95 and er 8 ~ 0.8/i _1 Mpc today for /3 = 0. In 
the left panel of Fig. 2 we show the evolution of energy 
densities, f2j(a) as a function of scale factor a, for an un- 
coupled ((3 — 0) and a coupled ((3 = 0.1) model. The 
most apparent feature here is that the epoch of matter- 
radiation equality z eq occurs earlier in the coupled model. 
This is due to the continuous flow of energy from the 



CDM to the quintessence field, i.e., $<p < and so if we 
evolve each model to the same cosmological parameters 
today, the energy density of CDM for the coupled 2EXP 
model will be higher at earlier times relative to the un- 
coupled one. In the right panel of Fig. 1 we show the 
redshift evolution of the Hubble parameter H{z) for a 
coupled and uncoupled 2EXP model, which we compare 
to ACDM. 

For p < 0, a shallow minimum is introduced to the 
effective potential, allowing <p to change sign during 
its evolution. This minimum is shallow enough that 
the field only undergoes a single oscillation, inducing 
damped oscillations in the equation of state as 
can be seen from the left panel of Fig. 1. Hence, the 
direction of energy exchange between the two coupled 
dark components is reversed during cosmic evolution. 
When normalised to common cosmological parameters, 
the net result is that for f3 < the CDM density scales 
less rapidly than the uncoupled case and the energy 
density of CDM will be lower at earlier times. As a 
consequence of the minimum in the effective potential 
the 2EXP model with f3 < also demonstrates some 
early dark energy behaviour. 

Albrecht Skordis (AS) - Albrecht and Skordis [42] pro- 
posed a phenomenologically motivated polynomial coef- 
ficient in their scaling potential, 

V^) = V Q [A+{<P-B) 2 ]e~ x \ (8) 

which introduces a local minimum to the potential, in 
which the quintessence field can get trapped. This model 
is able to satisfy a number of cosmological constraints 
and provide late time acceleration, whilst keeping all the 
parameters in the potential ~ C(M p i). We set A = 8.0, 
A = 0.01 and B = 33.960 throughout this study, val- 
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FIG. 2. The cosmological evolution of the energy density parameters Qi in the 2EXP (left panel) and AS (right panel) models: 
fi 7 (a) (solid), f2 c (o) (short-dashed), £lb(a.) (long-dashed) and Q^(a) (dot-dashed). The black (red) lines are for /3 = 0(0.1). 
Note that z cq occurs earlier (later) in the coupled 2EXP(AS) cosmology as a result of its higher (lower) CDM density at early 
times. 



ues chosen to provide fi<£ ~ 0.7 today when /3 — and 
Vb ~ O(Mpi). With these parameter values, the field ini- 
tially has a large value of V(<p) and rapidly approaches 
an attractor solution where mimics the equation of 
state of the dominant background energy density. Soon 
after the field enters the matter era, where w$ should 
go to zero, the field reaches the minimum of the poten- 
tial and undergoes small damped oscillations about that 
minimum (see Fig. 1). Due to damping from the Hub- 
ble expansion, these oscillations eventually stop and the 
field settles in the minimum, behaving as a cosmological 
constant with w<p = — 1. These late time oscillations of 
the quintessence field are further damped when the field 
is coupled to cold dark matter with f3 > 0; the effective 
potential will have a much sharper minimum, in which 
case the field remains in the attractor regime for a longer 
period, before suddenly getting trapped in the minimum 
without oscillating. For progressively stronger couplings, 
the oscillations in the field become more damped. For 
all values of /3 considered in this study, the amplitude of 
the oscillations in the field are small enough that they 
have a negligible effect on the evolution of the CDM den- 
sity. For the AS model, H(z) is identical to ACDM over 
the redshift range of interest, for all coupling strengths 
considered in this study. 

In the right panel of Fig. 2 we show the evolution of 
energy densities for an uncoupled (f3 — 0) and a coupled 
(j3 = 0.1) model. Since f3cb > (the direction of energy 
transfer is from the quintessence scalar field to CDM) 
the energy density of CDM will be lower at earlier times 
relative to the uncoupled one if we evolve each model 
to the same cosmological parameters today. Hence, the 
epoch of matter-radiation equality z eq occurs later in 
the coupled model relative to the uncoupled one. Notice 
also the non-negligible contribution from dark energy 
at early times. With the above parameter values A, B 



and A in the potential, it is not possible to realise a 
cosmology with f2^ ~ 0.7 today for j3 < —0.05 and so we 
only consider j3 > for the AS models. 

Supergravity (SUGRA) - As the vacuum expectation 
value of the quintessence field is of the order the Planck 
mass, a number of authors [43, 44] have argued that 
realistic quintessence models should receive supergrav- 
ity (SUGRA) corrections. Such a model was proposed 
in [43] , with a potential of the form 

V(4>) = ^V /2 , 0) 

where the evolution of <fr is insensitive to a wide range of 
initial conditions. Throughout this work we use A = 11 
to drive close to —1 today for j3 = resulting in 
Vb ~ 10 _8 M p i. Initially the value of the field is small 
in comparison to the Planck mass, and so the exponen- 
tial factor only plays a role at low redshifts, where its 
presence reinforces the domination of the potential en- 
ergy over the kinetic one, which pushes Wrj, toward — 1 
today. This is demonstrated in Fig. 1 where we show the 
evolution of for an uncoupled and coupled (/3 = 0.1) 
model. 

In Fig. 3 we depict the evolution of energy densities 
for an uncoupled ((3 = 0) and a coupled ((3 = 0.1) 
SUGRA model. Here, z cq occurs later in the coupled 
model relative to the uncoupled one, since the evolution 
of /3(f> > 0, i.e., the direction of energy transfer is from 
the quintessence scalar field to CDM. For f3 < the 
situation is reversed; j3<f> < and the energy density 
of CDM will be higher at earlier times relative to the 
uncoupled with z cq occurring earlier. The SUGRA 
models also have a non-negligible DE contribution at 
early times. In the right panel of Fig. 1 we show the 



5 




G cff -G N (l + 2/3 2 ), 
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FIG. 3. The cosmological evolution of the energy density 
parameters Qi for the SUGRA potential. Line types as in 
Fig. 2. Notice that z eq occurs later in the coupled cosmology 
as a result of lower CDM density at early times. 



redshift evolution of the Hubble parameter H(z) for a 
coupled and uncoupled SUGRA model. 

In summary, the direction of energy transfer, that is, 
the sign of the source term Q( c )o = —fip c 4>-, is dictated by 
the sign of j3 and shape of the effective potential in which 
the quintessence field evolves. In a coupled model where 
< (> 0), p c scales less (more) rapidly than (3 = 
and therefore if we normalise to the same cosmological 
parameters today, the energy density of CDM will be 
higher (lower) at earlier times. This normalisation to 
common cosmological parameters at redshift z = allows 
a direct comparison between different models. 



B. Linear perturbations 

In this section we discuss perturbation evolution in 
the coupled 2EXP, SUGRA and AS models, comparing 
them to uncoupled models and to ACDM. 

In the matter dominated epoch and on sub-horizon 
scales, whilst perturbations to the uncoupled baryonic 
component, 5^ = <5pb/Pb, obey (cf. Ref. [14]) 



-H5' h + -H 2 [n c 6 c + Q h S h ] , (10) 



perturbations to the coupled CDM component, S c = 
6p c /p c , follow 



(11) 



with 



where Gn is Newton's constant. Here, the prime denotes 
differentiation with respect to conformal time and W = 
a' /a is the conformal Hubble function. In the expressions 
above we have assumed that the scalar field is massless, 
which is a very good approximation for the models we 
consider here, since the length scales we are interested 
in are much smaller than the interaction range of the 
quintessence field. 

For /J / 0, the growth of cold dark matter pertur- 
bations are affected in a number of ways [19, 45, 46]. 
Firstly, the background expansion rate H{z) is modi- 
fied, which will lead to faster or slower clustering of mat- 
ter particles. Secondly the friction term of Eq. (11) is 
modified by the extra contribution, /3</>', implying that 
the rate of perturbation growth will then depend both 
on the coupling /3 and the background evolution of the 
scalar field. Secondly, the CDM particles feel an addi- 
tional 'fifth-force' (an effective Newton's constant G e ff) 
which is mediated by the scalar field. Finally, the cou- 
pling effectively rescales the CDM particle mass which 
the uncoupled baryonic component will feel through a 
changing gravitational potential generated by the CDM. 

In Fig. 4 we show the evolution of the ratio of the 
CDM density contrast in the coupled 2EXP and SUGRA 
models to ACDM, <5 c /<5cAcdm- These figures were ob- 
tained by numerically integrating Eqs. (10) and (11). For 
each model, the same initial conditions 8 l c b and 5 l c b were 
used and the integration was started at a high redshift 
(l + Zi = 10 9 ) to minimise the effect of the decaying mode 
on the evolution of the perturbations. 

For the 2EXP model with (3 > 0, the (3(j)' term con- 
tributes negatively to the total friction term of Eq. (11) 
which, along with the additional fifth force and higher 
CDM density at early times, leads to an enhancement of 
growth with respect to ACDM. For (3 = 0, the presence 
of the scalar field leads to a suppression of perturbation 
growth. We also find a suppression of growth in the un- 
coupled SUGRA and AS models. This is a well under- 
stood feature of uncoupled quintessence [32, 47, 48]; just 
as fluctuation growth is suppressed at late times with the 
onset of dark energy, the growth of linear perturbation 
modes are slowed at early times due to a non-negligible 
fraction of quintessence scalar field energy density. 

Despite the presence of the fifth force, CDM perturba- 
tion growth in the SUGRA and AS models is suppressed 
further for (3 > with respect to ACDM due to the 
lower CDM density at early times and since f34>' con- 
tributes positively to the total friction term. For (3 < 0, 
growth in the SUGRA model is enhanced compared to 
the uncoupled case and exceeds ACDM for progressively 
more negative couplings, due to the combined effect of 
the fifth force and since the f3(j)' term is anti-frictional. 
This behaviour is also true of the AS model. For (3 < 
in the 2EXP models, the growth of CDM perturbations 
has a more complicated dependence on the magnitude of 




FIG. 4. Evolution of the ratio of the CDM density contrast to ACDM, 5 c /S c acdm, for the coupled 2EXP (left panel) and 
SUGRA (right panel) models. These figures were obtained by numerically integrating Eqs. (10) and (11). Notice that the 
scales are different. 



/3, due to the reversed direction of energy exchange be- 
tween the CDM and quintessence field, induced by the 
shallow minimum in the effective potential. For weak, 
negative couplings /3 ~ —0.05, the net effect is a sup- 
pression of growth relative to ACDM and the uncoupled 
2EXP model. When the coupling strength is increased to 
j3 ~ —0.1, (3<j> < as the field rolls toward the minimum 
of the effective potential. This contributes negatively 
to the total friction term of Eq. (11), which combined 
with the fifth force acts to slightly enhance perturbation 
growth at early times (z ~ 100) with respect to ACDM 
and the uncoupled 2EXP model. Once the field turns 
over, /3(f) > in competition with the fifth force, with 
the net result that perturbation growth is suppressed at 
late times (z ~ 20) relative to ACDM and the uncoupled 
2EXP model, but enhanced relative to j3 — —0.05. For 
even stronger coupling strengths (3 ~ —0.2, the effect of 
the small frictional term at early times and the fifth force 
result in a enhancement of growth relative to ACDM. In 
practice, the impact of these individual effects on cosmic 
structure formation is not always obvious; one mecha- 
nism may act to enhance structure formation, whilst an- 
other may suppress it (cf. Ref. [46]). 

It is clear that changes in the background expansion 
history of the universe and modifications to the growth 
rate of the matter perturbations will affect the matter 
power spectrum. In Fig. 5 we plot the linear (total) 
matter power spectrum at z = 0, Po(fc), generated using 
CAMB. All the spectra have been normalised to CMB 
fluctuations. The left, middle and right panels show the 
power spectra for the 2EXP, AS and SUGRA models 
respectively for a variety of positive coupling strengths 
j3. For the 2EXP model, when the coupling strength is 
increased, small scale power is enhanced and a slight 
suppression of power on large scales is observed. This is 
because the epoch of matter-radiation equivalence (z cq ) 
occurs earlier when /3 > and so only the very small scale 



perturbations on scales k > k oq have time to enter the 
horizon and grow during radiation domination. Hence, 
the turnover in the power spectrum (which corresponds 
to the scale that entered the horizon when the universe 
became matter dominated) occurs on smaller scales and 
since the growth of perturbations for j3 > are en- 
hanced, small scale power is increased. For fluctuations 
on scales k < k cq during the matter dominated epoch, 
the enhancement or suppression takes place after the 
mode enters the horizon. Consequently, the amplitude 
of the linear power spectrum on 8ft. _1 Mpc _ scales (a$) 
is higher compared to ACDM, whilst the uncoupled 
2EXP model has a lower erg compared to ACDM when 
normalised to CMB fluctuations. Similar arguments 
may be applied to explain the suppression/enhancement 
of power seen on different scales for the AS and SUGRA 
cosmologies (middle and right panels of Fig. 5) . 



III. ABUNDANCE OF CDM HALOS 

While N-body simulations are a powerful tool for 
studying cosmic structure formation, they can be very 
time consuming. Therefore semi-analytic methods are 
a useful alternative, in particular for exploring a range 
of non-standard cosmological models involving a wide 
parameter space. The semi-analytical expressions used 
below have been shown to be in reasonable agreement 
with N-body simulations at calculating the abundance of 
dark matter halos over a wide range of masses, however 
we note that the expressions for the mass functions used 
below have been tested only against N-body simulations 
of ACDM cosmologies and not for coupled quintessence 
models. Fortunately, N-body simulations presented in 
[19] showed evidence that the Jenkins et al. mass func- 
tion [28] represents a good fit to the mass functions found 
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FIG. 5. The linear (total) matter power spectrum at z = 0, Po(k), generated using the CAMB computer software for the models 
under consideration. We show various coupling strengths and compare to the linear power spectrum of ACDM in the lower 
panels. Left panel: 2EXP model. Middle panel: AS model. Right panel: SUGRA model. Note that all spectra have been 
normalised to CMB fluctuations, which are on larger scales than included in the figure. 



numerically even in coupled quintessence models. 



A. Mass functions 

Observations indicate that individual galaxies and 
clusters of galaxies are embedded in extended halos of 
dark matter. Press and Schechter [29] were the first to 
provide a theoretical description for the abundance of 
these halos as a function of their mass using the assump- 
tion that the fraction of the volume of the universe that 
has collapsed into objects of mass M at a redshift z fol- 
lows a Gaussian distribution: 



J*col(M, z) 



'2ira(R, z) Js t (z) 



dSe -5*/2*HR,z) ^ (13) 



where 8- k {z) is the critical density which an object must 
exceed for collapse to occur at the redshift z of interest. 
The density field has been smoothed on a comoving scale, 
R, and <r 2 (R, z) is the mass variance: 



a 2 (R,z) = 



X*) 



2ir 2 



dkk 2 P {k)W 2 (k-R), (14) 



where Po(k) is the linear matter power spectrum today 
and W(k;R) is the Fourier transform of the window func- 
tion used to smooth the density field. The redshift de- 
pendence of the power spectrum has been extracted by 
multiplying by the linear growth factor of the perturba- 
tions of the total matter (p m = p c + pb) perturbations, 
normalised to unity today: 



D+(z) 



6 m (z = 0) ' 



(15) 



where S m = Sp m /p m . For a spherical top-hat filter of 
radius r in real space, containing a mass 



M (r, z) 



4 ft„(z) 

— 7T 1 

3 {l + zf 



(16) 



where p m (z) is the mean total matter density at redshift 
z, then 



W(k;R) 



_ /sin(fci?) cos(fci?)\ 
~ ! IkRY {kRy J 



(17) 



To obtain the comoving number density of objects of 
mass between M and M + AM at redshift z, we first 
differentiate Eq. (13) with respect to M which gives the 
fraction of the volume in the universe that has collapsed 
into objects with mass between M and M + dM. Mul- 
tiplying this by p m (z)(l + z)~ 3 /Al gives the comoving 
number density of collapsed objects with mass between 
M and M + dM at redshift z: 



dn(M, z) _ 1 Pm(z) dT co \ 
dM ~ M(l + z) 3 dM 

Carrying out the differentiation gives 



(18) 



dn(M, z) 
dM 

x exp 



2 1 p m (z) 8*(z) 
7T M (1 + z) 3 a(R, z) 



5 2 Xz) 
2a 2 {R,z) 



1 da{R,z) 
~a(R,z) dM 



(19) 



and is referred to as the differential mass function. Writ- 
ing the term in the square brackets as the logarithmic 
derivative d lner^/dA/, it is common practice to re- 
write Eq. (19) as 
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(l + z) 3 dn(M,z) 

/ K z) = M _ — — 

p m (,z) diner 1 



(20) 



which defines the mass function. Hence, (dropping the 
arguments of 5±(z) and a{R 1 z) for compactness) the PS 
mass function is 



/(<t;PS) = a/ 2 ^ exp 

V 7T (T 



2ct 2 



(21) 



and so has an implicit dependence on redshift, through 
the parameter 5+(z). The PS mass function is normalised 
so that 



/(cr; PS) diner" 1 = 1, 



(22) 



which implies that all of the matter is within haloes of 
some mass. It has been found to over-predict (under- 
predict) the number of low(high) mass haloes at the cur- 
rent epoch [28, 49]. Shcth and Tormen [30] (hereafter 
ST) found that this discrepancy can be substantially re- 
duced by considering ellipsoidal, rather than spherical 
collapse. Their mass function contains three empirically 
determined parameters A, a and p: 



/fa ST) 



1 



a5l 



<5* 



exp 



2cr 2 



(23) 

We use A = 0.322, a = 0.707 and p = 0.3 [30], for all the 
quintessence models considered. The functional form of 
/ in Eq. (20) has also been directly fitted to the mass 
function obtained from numerical N-body simulations as 
a function of a only. For example, Jenkins et al. [28] 
found 

/(cr; Jenkins) = 0.315 exp (-| In cr" 1 + 0.61| 3 8 ) , (24) 

valid over the range —1.2 < In cr -1 < 1.05. Note the 
range over which the fit function is valid. For standard 
cosmologies, this corresponds approximately to the mass 
range 10 11 to 10 16 M Q at z — 0. We find that for some 
combinations of quintessence potential and coupling, 
the mass/redshift range over which this fit function is 
valid is dramatically reduced. Since Eq. (24) depends 
on the variance a only, it is manifestly independent 
of cosmology-specific parameters such as redshift and 
critical density for collapse. This is what is commonly 
understood as the 'universality of the mass function', 
(see Refs. [50] and [52] for recent discussions of this 
topic). Mass function universality suggests that semi- 
analytic forms such as ST and fit functions like Jenkins 
should remain valid in coupled quintessence cosmolo- 
gies. Indeed, this has been verified by a number of 



authors [19, 20] by using N-body simulations. Another 
mass function, similar in form to ST has been introduced 
by Tinker et al. [51]. It has been calibrated over a wide 
mass and redshift range using N-body simulations, 
but is however, still only valid over a limited range in 
In cr -1 . Hence, for our non-standard CQ models, the 
Tinker mass function offers no real advantages over 
ST, since to probe high (z > 1) rcdshifts would require 
extrapolating the Tinker fit-function beyond its range 
of validity. Recently, a path-integral method that allows 
for an analytical computation of the mass function for 
more generic filters than the one described above, has 
been used to provide a robust theoretical estimation 
of the halo mass function [52]. Based on the Path 
Integral approach applied recently to the Excursion Set 
formalism of Ref . [53] , these theoretical predictions show 
a remarkable agreement with the numerical results of 
Tinker et al. with deviations of no greater than 5% over 
the range of masses probed by the simulations. This 
suggests that the Excursion Set formalism, combined 
with a realistic model of the dark matter halo collapse 
conditions can provide a universal formulation of the 
halo mass function. Although, not the thrust of this 
piece of work, it would certainly be worth following up 
this possibility for the case of coupled quintessence. 

A cumulative number count Af(> M, z, Az), above a 
given mass threshold M — M th , in a redshift slice Az is 
obtained from dn(M, z)/dM: 



f z+Az dV 
Af(> M, z, Az) = AO J dz—n(>M,z 



n(> M,z) 



I 



- dM MM,z) 

M th 



dM 



(25) 



(26) 



where dV/dz is the comoving volume element, which in 
a spatially flat universe takes the form 



dz= AQr{z) dz 



(27) 



AO is the solid angle and r(z) denotes the comoving ra- 
dial distance out to redshift z: 



r(z) = 



dx 

H Jo W) 



(28) 



where E(z) = H{z)/Ho- This geometrical factor dV/dz 
is required since the redshift evolution of a physical 'vol- 
ume' in space is model dependent, depending (among 
other things) on the coupling strength and form of the 
potential. For numerical computation, the upper limit of 
integration in Eq. (26) is replaced by some finite mass 
value M max . 
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The final ingredient required is a method for calculat- 
ing the critical density contrast for collapse, 5*(z). In 
a matter dominated universe, this parameter is a con- 
stant 5* = 1.686, whilst for ACDM it becomes a weak 
function of cosmology, asymptoting to <5* = 1.686 at 
high redshift [54, 55]. In cosmologies with an early dark 
energy component however, it has been established (see 
Refs. [56, 58] for example) that the value of 5+ can have a 
small, but non-negligible effect in the high mass/redshift 
end of the mass function, where the cluster abundance 
is smallest. As discussed in e.g. Refs. [57] and [58], if 
DE were not to remain homogeneous over the relevant 
length scales, or if, as is the interest here, DE is coupled 
to dark matter, then the DE density within a pertur- 
bation may evolve separately from the background uni- 
verse. We should then expect significant deviations from 
5* = 1.686. Determining S+(z) in CQ cosmologies has re- 
cently been the subject of some debate [31]. In the next 
section, we present the spherical collapse model applied 
to CQ and the model proposed in Ref. [31]. Both meth- 
ods are based on assumptions: In Ref. [31] for example, 
one neglects the contribution of the scalar field pertur- 
bation in the Poisson equation, whilst in the spherical 
collapse model one assumes that there is no loss in scalar 
field energy density from inside the collapsing region. To 
get a feeling for how the lack of understanding in mod- 
elling the non-linear behaviour of the scalar field man- 
ifests itself in the halo mass function, we subsequently 
compare these two methods for determining 5+(z). 



-2H5 m + lp m 6 m (l + S m ) + 

2 S 1 + o n 

1 



+ ^ [M 1 + 3^0) - + 3^)] (1 + <5 m ) 
+ (2H5'<t> + ty) (1 + S n 



2 

+13 



2 
3 

1 



+p 2 -6f(l + 5 m ). 



(29) 



The evolution equation for the scalar field inside the over- 
density reads: 



3 \H+^ 



1 S n 



31 + S n 



dV 



= -/3p m (l + S n 



B. Critical density for collapse, <5*(z) 



(30) 

From Eq. (29) we see that the original non-linear equa- 
tion for the density contrast is modified by the presence of 
the coupling terms and by a term that originates from the 
assumption that the quintessence field inside the over- 
dense region evolves separately from the background. 

Since we are interested in calculating the extrapolated 
linear overdensity of a matter perturbation collapsing at 
any given redshift, i.e., 5±(z) = 5 m x(z = z c ), we solve 
numerically Eqs. (29) and (30) with 5 % m = and search 
for the value of the initial overdensity S l m for which the 
collapse occurs at redshift z c . Formally, this is when 
Sm{z c ) — > oo, i.e., the sphere collapses toward a central 
singularity. Of course, in reality, dissipative processes 
intervene well before such a singularity is reached and 
the system will virialise. Having determined the initial 
overdensity for the redshift of interest, we use it as an 
initial condition to solve the linearised form of Eq. (29) 
which defines the linearly extrapolated density contrast 

<5m,L, 



In this section, we denote background quantities with 
a bar, we assume matter domination and we neglect 
baryons. The evolution of a density perturbation during 
the non-linear phase of gravitational collapse is most sim- 
ply approximated by the spherical collapse model [59, 60]. 
The model was first applied to the standard CDM sce- 
nario and later to ACDM [54]. Since then, the spherical 
collapse model has been applied to quintessence scenar- 
ios [57, 61, 62] and more recently to interacting DE-CDM 
models [63]. In an alternative method, Ref. [21] made 
predictions for 6+(z) in CQ models, explicitly accounting 
for an uncoupled baryonic component by considering the 
evolution of spherical concentric shells. 

All background quantities (H, (f>, p m ) follow the stan- 
dard evolution equations for coupled quintessence (cf. 
Ref. [14]). Defining SJ> as the perturbation to the 
quintessence field, (f> = <f> + 6<f>, the non-linear evolution 
of the density contrast is governed by [63] 



+ \ [(! + 3u^)p^ + 3p$6wj] 

+P(2H6<f> + 8<j)) , (31) 

where we have defined p = p^(l + 6$), 6$ = Sp^/p^. 
Furthermore, we have 

= a . ( 32 ) 

and 

dV 

5 P j = 4>6<l>+-^6<f>. (33) 

In this linear regime the quintessence field inside the over- 
density obeys 
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FIG. 6. The redshift evolution of the linearly extrapolated density contrast at the redshift of collapse <5*. We show various 
coupling strengths (3 > and compare to the ACDM value (thick solid back line). Left panel: 2EXP model. Middle panel: 
AS model. Right panel: SUGRA model. Notice that the scales are different. In all panels, thick(thin) lines correspond to the 
Newtonian limit (spherical collapse model). Within the framework of the PS and ST formalism, cosmologies with lower values 
of 8^ form structure earlier. 



3[H 



-PPm(l 



(34) 

Eqs. (30-34) were first derived in Ref. [63] for the case 
of a variable coupling /3(4>). In arriving at this descrip- 
tion for the non-linear evolution of a spherical matter 
perturbation, approximations have been made. Firstly, 
it is assumed that there is no loss in scalar field en- 
ergy density from inside the collapsing region [57, 63] 
- the quintessence field is assumed to cluster with the 
CDM and whole system (CDM and scalar field) virial- 
izes. It has been argued [31] that since the standard 
spherical collapse model is intrinsically based on gravi- 
tational attraction only, it incorrectly describes the evo- 
lution of overdensities into the non-linear regime when 
applied to CQ cosmologies. Indeed, when a dark-dark 
coupling is introduced, CDM will experience an addi- 
tional attractive force which is not explicitly included in 
the spherical collapse formulation. By considering the 
Newtonian limit of the relativistic non-linear perturba- 
tion equations, Ref. [31] was able to show that the evolu- 
tion of the non-linear and linear density contrasts obey 



5 m = -(2H + /3(/ ) )d n 
and 



4 5 2 



3 1 



+ ^p m 5 m (l + 2/3 2 ), (35) 



1 



<5 m , L = -(2H + W)<5 m , L + -p m 6 m ±(l + 2/3") . (36) 

respectively. We refer the reader to Ref. [31] where a 
detailed comparison to the spherical collapse model is 
presented. Eq. (36) is simply the relativistic Eq. (11) in 
standard cosmic time. We immediately notice (as the 
Newtonian limit demands) that Eqs. (35) and (36) con- 
tain no terms involving time derivatives of 8<j). Also, there 



is no term analogous to the fourth and third term on the 
RHS of Eqs. (29) and (31), since to satisfy the Newtonian 
limit the quintessence field should not contribute to the 
gravitational potential. Notice that this model reduces 
to the standard spherical collapse description in the limit 
of a cosmological constant with zero coupling. As we will 
now see, this approach can lead to significantly different 
predictions for S+(z) when compared to the spherical col- 
lapse model. 

In Fig. 6 we depict the redshift evolution of the 
parameter 5±(z), where the thick lines represent the 
Newtonian (NWT) prediction and the thin lines the 
corresponding spherical collapse (SPH) prediction. We 
depict a range of coupling strengths j3 > 0. Notice 
the SPH model prediction of oscillatory behaviour of 
<5* (z) for the uncoupled AS model. This is due to late 
time oscillations of the field about the minimum of the 
potential and are progressively damped for stronger 
couplings (see Section II A for an explanation). These 
oscillations are not present in the NWT prediction since 
the quintessence field is assumed not to cluster. In 
general, for the 2EXP potential (for which perturbation 
growth is enhanced relative to ACDM for j3 > 0), the 
NWT model predicts a higher 6+(z) with higher j3, whilst 
the SPH model predicts lower values of 8*(z) with higher 
(3. For the SUGRA and AS potentials (where perturba- 
tion growth is suppressed relative to ACDM for j3 > 0), 
the NWT model predicts a lower 6+(z) with higher /?, 
whilst the SPH model predicts higher values of 6+(z) 
with higher (3. For /3 < 0, this behaviour is reversed. 
Within the framework of the PS and ST formalism, 
cosmologies with lower values of cv form structure earlier. 

In summary, the assumptions made when following 
the dark energy scalar field into the non-linear regime, 
i.e., whether the field is to be treated separately inside 
and outside an overdense region, greatly effects the crit- 
ical density at the redshift of collapse. Without proper 
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knowledge of how to treat the non-linear evolution of 
the field, we compare halo abundance computed using 
the NWT and SPH approaches. The differences between 
the SPH and NWT models for some CQ cosmologies are 
large enough to become manifest in the high-z, high- 
mass tail of the mass function due to the exponential 
dependence on It is important however to put these 
differences into context with the other uncertainties that 
go into modelling the mass function. As we will dis- 
cuss in section IV B, the differences between the NWT 
and SPH models are secondary compared to the uncer- 
tainties associated with the use of the extremal values of 
/3, obtained by fitting the quintessence models to data. 
When one then considers the limited accuracy of the fit- 
ting functions f(a, z) used in this work, the difference 
induced in the mass function by using either the SPH or 
NWT models becomes a secondary concern. 



IV. RESULTS 

We will now present results for the halo abundance 
for the coupled quintessence models introduced in 
Sec. II A using the mass functions described in Sec. III. 
To facilitate an easier comparison between models, we 
begin by normalising each model to the ACDM WMAP7 
Maximum Likelihood (ML) cosmological parameter 
values [2]. As the resulting cosmologies are unlikely to 
be consistent with all observational data, as pointed 
out in Ref. [32], we subsequently perform a global fit 
using the CosmoMC package [33], a Monte Carlo Markov 
Chain (MCMC) code, to find the best fit values for the 
cosmological parameters such that each model satisfies 
the CMB+BAO+SNla+ifo datasets. 



A. WMAP ACDM cosmological parameters 

In this subsection we vary the normalisation of the 
potential and the CDM density, so that all cosmologies 
share the same ACDM WMAP7 cosmological parameter 
values [2], given in Table I of Appendix A. For all un- 
coupled cosmologies (i.e. f3 = 0), we find a lower abun- 
dance of high mass haloes compared to ACDM, when 
normalised to CMB fluctuations. This is in agreement 
with previous studies (c.f. [64] and references therein). 

In Fig. 7 we show the theoretically expected abun- 
dance of CDM haloes, J\f(> M, z, Az), in the mass range 
10 n M Q /i- 3 < M < 10 15 /j- 3 M© and in a redshift vol- 
ume z = 0.01—0.4, using the Jenkins et al. mass function. 
This mass range extends from typical galaxy size halos, 
up to galaxy cluster masses. Below each figure we plot 
the ratio of the quintessence cosmology to ACDM to facil- 
itate an easier comparison. Starting with the AS model, 
the Jenkins et al. mass function predicts a reduction in 
cluster abundance in the high mass tail. Numbers of the 
most massive clusters are further reduced when a cou- 



pling j3 > is introduced. The same trend is true of the 
SUGRA models, however suppression in the high mass 
tail is more severe, with differences compared to ACDM 
reaching almost two orders of magnitude for haloes of 
mass M > 1O 14 /i~ 3 M for /3 = 0.1. Here, the sign of 
the coupling parameter and the shape of the effective 
potential in which the field evolves dictates that there is 
a continual flow of energy from the quintessence field to 
the CDM, i.e., <3( c )o < 0. Considering the arguments of 
Sections II A and II B, the effect of this direction of en- 
ergy exchange is to suppress structure formation relative 
to ACDM, regardless of the fifth force mediated by the 
quintessence field. 

As one moves from ACDM to the uncoupled 2EXP 
model, we expect to see a reduction in the number of 
high mass (M > 10 14 ft.~ 3 M Q ) clusters. When a weak 
positive coupling (/3 = 0.05) is introduced, ACDM abun- 
dances are recovered, and as the coupling strength is in- 
creased we expect to see an excess of high mass clusters, 
since Q( c )o > 0. In the SUGRA model, for increasingly 
more negative couplings, Q( c ) > and ACDM abun- 
dances are recovered and then surpassed. In the 2EXP 
model, for weak negative couplings j3 ~ —0.05, the num- 
ber of high mass clusters is suppressed relative to (3 = 0. 
For stronger, more negative couplings, high mass clus- 
ter numbers in the uncoupled model are recovered and 
with (3 <^ —1.2, ACDM numbers are surpassed. These 
effects are demonstrated in the bottom-right panel of 
Fig. 7 where we compare the three potentials for selected 
positive and negative coupling values. As has been no- 
ticed previously [21], the differences between the uncou- 
pled quintessence models and ACDM can be partially or 
even totally reduced, (i.e. they become degenerate), by 
introducing a weak DE-CDM coupling. 

In the upper panels of Fig. 8 we compare the halo 
abundances calculated using the PS and ST mass func- 
tions to the Jenkins et al. mass function, for the 2EXP 
and SUGRA models with f3 = 0.1. We use both the NWT 
and SPH methods for calculating 5*(z). Starting with the 
2EXP cosmology, ST-NWT agrees well with the Jenk- 
ins prediction, whilst ST-SPH predicts a slight excess 
of high mass haloes compared to both Jenkins and ST- 
NWT. The PS mass function predicts a slightly greater 
abundance of low mass clusters compared to Jenkins. At 
the high mass end, the PS-SPH model lies above the 
Jenkins prediction, whilst the PS-NWT predicts a reduc- 
tion in halo abundance. The same general trend between 
the PS, ST and Jenkins et al. mass functions is true of 
the SUGRA cosmology, however this time the spherical 
collapse model predicts a lower number of haloes com- 
pared to the Newtonian model. These differences are 
because in the 2EXP cosmology, the NWT model pre- 
dicts a higher 5 it (z) with higher /3, whilst the SPH model 
predicts lower values of 5* (2) with higher f3. Within the 
framework of the PS and ST formalism, cosmologies with 
lower values of 5* form structure earlier, and hence the 
SPH model predicts a greater abundance of haloes. The 
same argument applies to the SUGRA cosmology, where 
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FIG. 7. The Jenkins et al. mass function prediction for the mass distribution of CDM haloes with 10 11 /i _3 Mq < M < 
10 15 /i _3 Mq, within a redshift survey volume z — 0.01 — 0.4. Top left panel: AS model. Top right panel: SUGRA model. 
Bottom left panel: 2EXP model. Bottom right panel: A comparison between the three quintessence models for a range of 
coupling values. Note that the Jenkins et al. mass function does not reach the highest cluster masses, since it is only valid over 
a limited range of a. 



for /3 = 0.1, the NWT model predicts lower values of 
8±(z) compared to the SPH model. For the AS cosmol- 
ogy, we found that the differences between the NWT and 
SPH model predictions for S i ,(z) were negligible when in- 
corporated into the PS and ST mass functions. 

Focusing on the high mass tail of the mass function, 
since this is where deviations from ACDM are most 



prominent, we found that for stronger coupling strengths 
the 2EXP cosmology has a higher abundance of massive 
objects at high redshift compared to ACDM (for a given 
mass function), whilst deviations from ACDM remain 
much less than one order of magnitude at low z < 0.5 
rcdshifts. The uncoupled and coupled SUGRA and 
AS cosmologies have fewer massive objects at high and 
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FIG. 8. A comparison between the cluster number counts, N[> M,z,Az), for the different mass functions for lO n /i" 3 M < 
M < W 15 h~ 3 M Q , within a redshift survey volume z = 0.01 - 0.4. Top left panel: 2EXP model, j3 = 0.1. Top right panel: 
SUGRA model, j3 — 0.1. Below each figure the PS and ST predictions are compared to the Jenkins et al. mass function 
by taking the ratio. Note that the Jenkins et al. mass function does not reach the highest cluster masses, since it is only 
valid over a limited range of a. Bottom left panel: The redshift evolution of the cluster number count, J\f(> M,z,Az), for 
10 14 /i" 3 M Q < M < 10 15 /i" 3 M Q in the 2EXP model for (3 = 0, with redshift bins of width Az = 0.1. Bottom right panel: 
2EXP, /3 — 0.1. Below each panel, we plot the ratio of halo abundance to that found using the PS-NWT approach. 



low redshift and this reduction is enhanced for stronger 
couplings. In the lower panels of Fig. 8 we show the 
redshift evolution of the cluster number count for objects 
of mass 10 14 h~ 3 M Q < M < I0 15 h~ 3 M in the 2EXP 
model. For the uncoupled model (left panel), the NWT 



and SPH predictions for 8* agree well over the redshift 
range of interest, with PS theory underestimating the 
number of high redshift objects compared to the ST 
mass function. For the /3 = 0.1 model (right panel) 
the different redshift evolution of 5+(z) predicted by 
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FIG. 9. Left panel: The number counts of CDM haloes with lO 11 ^ 3 M Q < M < I0 15 h~^ 3 M s , within a redshift survey volume 
z — 0.01 — 0.4 found using the Jenkins et al. mass function. Right panel: The redshift evolution of the cluster number count, 
JV{> M,z,Az), for 10 14 /iX 3 Mq < M < W^Ii^Mq using the ST-NWT mass function. h A is the value of h in the ACDM 
cosmology. We use bins of redshift width Az = 0.1. 



the NWT and SPH models is clearly manifest at high 
rcdshifts. Similar, non-negligible differences occur 
between the NWT and SPH models in the SUGRA 
cosmologies. 



B. Maximum Likelihood cosmological parameters 

In this subsection we present mass function results for 
the 2EXP, SUGRA and AS models using the ML cosmo- 
logical parameter values obtained from a MCMC global 
fit. 

We take the set of base parameters (see Appendix A 
for an explanation of these parameters) 

ti = {VL h h 2 , Q c h 2 , H ,t, log [10 10 A|(fc )], n s , /?} , 

(37) 

and use a modified version of the CosmoMC package [33] 
to sample from the joint posterior distribution of the pa- 
rameters, 



£(a;|fl)P(fl) 
fd0£(x\0)P('d) 



(38) 



where C(x\-d) is the likelihood of the data x given the 
model parameters i? and "P(#) is the prior probability 
density. 



We use the CMB+BAO+SNla+-ff datasets (see Ap- 
pendix A) and apply a Gelman-Rubin [65] convergence 
criterion of R — 1 < 0.03 across eight chains for each 
model. We assume purely adiabatic initial conditions 
with no tensor contribution. We do not vary the pa- 
rameters X n in the quintessence potentials; they remain 
fixed with the values given in Sec. II A. We take this 
approach as our goal is to investigate the observability of 
couplings of order (3 ~ 0.1. Finding a best fit including 
varying the potential parameters is therefore not neces- 
sary for this particular approach although of course there 
is nothing to stop them being included in general. This 
approach also allows us to compare our results with those 
from Sec. IV A, and assess the effect of (erroneously) us- 
ing the ACDM cosmological parameters. For the 2EXP 
and SUGRA models we consider couplings in the range 
/3 € [—0.4, 0.4] i.e. including negative and zero cou- 
plings, whilst for the AS model, j3 G [0.0, 0.4], since it is 
not possible to realise a cosmology with ~ 0.7 today 
for f3 < —0.05. We assume flat priors (fix = 0) on all 
base parameters. In Table I of Appendix A we list the 
ML values of this parameter set and the derived param- 
eters for ACDM, 2EXP, SUGRA and AS. Also shown in 
brackets are the extremal values of the parameters in the 
full n-dimensional 68% confidence region. 

In Fig. 9, using the parameter values in Table I of Ap- 
pendix A, we show halo abundance as a function of mass 
(left panel) and as a function of redshift (right panel). 
The positive value of the coupling parameter (f3 = 0.033) 
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for the 2EXP-ML cosmology acts to slightly enhance 
the formation of the most massive haloes, which com- 
bined with the fractionally slower Hubble expansion rate 
(Hq = 69.86 km s -1 Mpc -1 ) and higher matter density 
(Q m = 0.273) mean that more massive objects have time 
to form during the matter era in the 2EXP-ML model 
relative to ACDM. These differences do not exceed a fac- 
tor of ~ 2 across the entire mass and redshift range of 
interest. Although the fractionally slower Hubble expan- 
sion rate (H = 69.5 km s" 1 Mpc" 1 ) for the SUGRA-ML 
cosmology will act to enhance structure formation, the 
lower matter density fraction (51 m = 0.221), negative 
coupling parameter (j3 — —0.078) and suppressed lin- 
ear power spectrum amplitude across all relevant scales 
(see Fig. 10) mean that less massive objects have time 
to form during the matter era in the SUGRA-ML model 
relative to ACDM, with abundances further decreasing 
with redshift, reaching one order of magnitude difference 
by z — 1.5. The slightly higher matter density of the 
AS-ML model (fi m = 0.273) compared to ACDM slightly 
raises the halo abundance, however the suppressed power 
spectrum amplitude (see Fig. 10) still means that num- 
bers of the most massive high-z clusters still fall short of 
ACDM predictions. Comparing these results with those 
from Section IV A, (where each model shared the same 
cosmological parameter values as ACDM) we see that the 
departures from ACDM halo abundance are not as great 
in the ML case. As the CQ models are by construction 
such that they are now consistent with a wider range of 
observational data, this is not unexpected. 
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FIG. 10. The linear total matter (p c + pb) power spectrum at 
z = 0, Po(k), for the three quintessence cosmologies with ML 
cosmological parameter values. Note that all spectra have 
been normalised to CMB fluctuations, which are on larger 
scales than included in the figure. 



V. CONCLUSIONS 



The differences in the threshold of collapse, 8±, ob- 
tained from the spherical collapse model and from the 
Newtonian limit of the non-linear perturbation equa- 
tions for the AS-ML model had a negligible effect on 
the halo counts when incorporated into the mass func- 
tion. We found non-negligible differences between the 
two approaches when used for the SUGRA and 2EXP- 
ML models, where the lower value of 6+ obtained from the 
spherical collapse model leads to higher halo abundance 
relative to the Newtonian prediction. As both the spher- 
ical collapse model and the model of Ref. [31] are based 
on approximations, we can, to some extent, view these 
differences as another source of uncertainty when making 
predictions for galaxy and cluster number counts. This 
uncertainly must be put into context however, alongside 
the other uncertainties, such as parameter constraints on 
j3 and the accuracy of the semi-analytic form of the mass 
function. For example, for all coupled quintessence mod- 
els, we find the differences induced in the mass function 
by taking the la n-dimensional extremal values of the 
coupling parameter f3 given in Table I of Appendix A, to 
be larger than the differences induced in the mass func- 
tion by using either the NWT or SPH models. When one 
then considers the limited accuracy of the fitting func- 
tions, f(a,z), used in this work, the modelling of the 
halo collapse threshold becomes a secondary concern. 



We have studied the abundance of CDM haloes in three 
different coupled quintessence cosmologies, using three 
common analytic forms for the halo mass function. For 
all uncoupled cosmologies, we found a lower abundance of 
high mass haloes compared to ACDM, when normalised 
to CMB fluctuations, which is in agreement with previous 
studies (cf. [64] and references therein). 

When the quintessence field is coupled to CDM 
through a coupling given by Eq. (5) the direction of en- 
ergy transfer is determined by the sign of the coupling 
parameter and on the effective potential in which the 
field evolves. For the models considered in this work, if 
there is a continual flow of energy from the quintessence 
scalar field to the CDM component, then the predicted 
number of CDM haloes can only lie below that of the un- 
coupled model (and therefore lower than ACDM) when 
each model shares the same cosmological parameters to- 
day. If the transfer of energy is from the CDM component 
to the quintessence field, then the predicted number of 
CDM haloes can be above or below ACDM depending on 
the strength of the coupling. Ensuring that the models 
are consistent with observational data by conducting a 
MCMC maximum likelihood analysis, the SUGRA-ML 
and AS-ML models predict reduced numbers of massive 
objects at high redshift relative to ACDM, whilst the 
2EXP-ML model yields an excess of the most massive 
haloes, increasing with redshift. Since the cosmological 
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parameter values obtained from the maximum likelihood 
analysis are different for each model, the predicted halo 
number counts for each model can differ from the counts 
obtained when each model shares the same cosmological 
parameter values today, for the same coupling value (3. 
This reinforces the point made in [32] that it is impor- 
tant to work with observationally consistent cosmologi- 
cal models when making predictions for cluster number 
counts, since the CDM halo mass function is exponen- 
tially sensitive to parameters such as fi m . 

It is intriguing that a coupled quintessence model can 
provide a cosmology that is consistent with observational 
data and produce an excess of high mass haloes. In light 
of the recent detections of massive high-z galaxy clus- 
ters [24-26] , a weakly coupled 2EXP quintessence model 
may help to alleviate tensions between these observations 
and the standard ACDM paradigm, for which these ob- 
jects should be extremely rare. This also suggests that 
CQ models where the coupling function is time depen- 
dent, 0(t), may serve as a viable alternative to ACDM 
(see also [20]). 

We have also shown that the assumptions made when 
modelling halo collapse in coupled dark energy cosmolo- 
gies (NWT or SPH), reflects strongly in the high-z, 
high-mass tail of the mass function due to the exponen- 
tial dependence on the threshold of collapse parameter 
S 2 . When this uncertainty is put into context however, 
alongside the other uncertainties such as parameter con- 
straints on (3 and the accuracy of the semi-analytic form 
of the mass function, the modelling of the halo collapse 
conditions becomes a secondary concern. That said, re- 
cent advances in mass function theory [52, 53], have al- 
lowed for a more robust calculation of the mass function. 
This formalism however still relies on a model for the 
halo collapse conditions. If it were applied to coupled 
quintessence, the approximations made when calculating 
the threshold of collapse, 5±(z), would become an impor- 
tant contributor to the total error budget. It will become 
important to be able to model halo collapse correctly in 
order to make best use of this new approach. With the 
arrival of greatly improved high redshift cluster surveys 
such as the Dark Energy Survey (DES) , Planck and even- 
tually LSST, the mass distribution of the most massive 
objects in the Universe will act as a powerful tool which 
may be used to constrain coupled quintessence models. 
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Appendix A: Maximum Likelihood parameter values 

Table I gives the Maximum Likelihood (ML) cosmo- 
logical parameter values for each coupled quintessence 
cosmology, obtained by performing a global fit using 
the CosmoMC package [33], a Monte Carlo Markov Chain 
(MCMC) code. We employed a Gelman-Rubin [65] con- 
vergence criterion of R — 1 < 0.03 across eight chains for 
each model. We assumed flat priors (CIk = 0) on all base 
parameters: 



= {n h h 2 , n c h 2 , H , t, log[10 10 A^(fc )], n s , 0) , 

(Al) 

and we assumed purely adiabatic initial conditions with 
no tensor contribution. Here, fl^h 2 and fi c /i 2 are the 
physical baryon and cold dark matter densities respec- 
tively, t is the reionization optical depth, A^(ko) is the 
primordial superhorizon power in the curvature pertur- 
bation on scales ko Mpc~ , n s is the scale spectral in- 
dex and f3 is the DE-CDM coupling parameter. For 
the 2EXP and SUGRA models, we consider couplings 
in the range /3 G [—0.4, 0.4], whilst for the AS model, 
j3 € [0.0, 0.4], since it is not possible to realise a cosmol- 
ogy with fl^ <~ 0.7 today for f3 <^ —0.05. The param- 
eters A„ in the quintessence potentials were kept fixed 
with the values given in Sec. II A. 

The observational constraints we use to evaluate the 
ML parameter values for each coupled quintessence cos- 
mology are local distance measures of the Hubble con- 
stant (Ho) and Type la Supernovae (SNIa) and preci- 
sion measurements of the Cosmic Microwave Background 
(CMB) and the Baryon Acoustic Oscillations (BAO). 
For the CMB, we use the 7-year data release from the 
WMAP satellite (WMAP7) [68] and utilise the likelihood 
code available at the LAMBDA website [69]. We com- 
pute the CMB angular power spectra using a modified 
CAMB code, described in section II A. The Hubble con- 
stant measurements are taken from the SHOES team [70] , 
based on geometric distance measurements of SNIa at 
0.023 < z < 0.1. We use the BAO constraints found in 
Ref. [71] and finally, we use the SNIa sample of Ref. [72] . 
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Class 


Parameter 


ACDM WMAP7-ML 


2EXP-ML (68%-full) 


SUGRA-ML (68%-full) 


AS-ML (68%-full) 


Primary 


won b h 2 

n c h 2 

H (kms^Mpc" 1 ) 
r 

log[10 10 A^(fc„)]f 

n s 
P 


2.227 
0.111 
71.4 
0.086 
3.170 
0.969 


2.244 (2.124 - 2.392) 
0.111 (0.101 - 0.121) 
69.86 (65.86 - 74.18) 
0.089 (0.052 - 0.131) 

3.08 (2.99 - 3.18) 
0.967 (0.941 - 1.01) 
0.033 (-0.169 - 0.096) 


2.329 (2.190 - 2.482) 
0.083 (0.076 - 0.910) 
69.50 (65.40 - 74.35) 
0.091 (0.055 - 0.140) 
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-0.078 (-0.123 - -0.044) 


2.281 (2.163 - 2.442) 
0.116 (0.108 - 0.127) 

71.4 (67.20 - 74.92) 
0.093 (0.055 - 0.134) 

3.07 (2.98 - 3.19) 
0.977 (0.946 - 1.010) 

0.0048 (0.0 - 0.046) 


Derived 


t (Gyr) 
O m 

z re 


0.738 
13.71 
0.262 
10.3 


0.727 (0.675 - 0.768) 
13.74 (13.37 - 14.09) 
0.273 (0.231 - 0.325) 
10.7 (7.4 - 13.9) 


0.779 (0.740 - 0.815) 
13.69 (13.24 - 14.0) 
0.221 (0.184 - 0.259) 
10.5 (7.4 - 13.7) 


0.727 (0.677 - 0.759) 
13.52 (13.22 - 13.89) 
0.273 (0.240 - 0.323) 
11.0 (7.7 - 14.2) 



TABLE I. Maximum Likelihood values for ACDM and the three coupled quintessence models. For the coupled quintessence 
we also show in brackets the extremal values of the parameters in the full n-dimensional 68% confidence region. ^) A 2 z (k) = 
k 3 Vn(k)/{2n 2 ). For ACDM WMAP7-ML, k = 0.002 Mpc" 1 . For 2EXP-ML, SUGRA-ML and AS-ML, k = 0.05 Mpc' 1 . 
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